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I. INTRODUCTION 



The physics of quantum phase transitions in one- 
dimensional systems has attracted much attention in re- 
cent years because not only of its speciality — the con- 
formal symmetries and universalities in (l+l)-critical 
modelsji*2iii but also of its generality — the quantum criti- 
cal phenomena in higher dimensions^ A general question 
is posed as follows: If a critical system defined in the 
ultraviolet limit is deformed by more than one relevant 
perturbations, what is the fate of the system? Can the 
criticality be reached again in the infrared limit because 
of the "cancellation" effect among these perturbations, 
thus a quantum phase transition takes place? A typical 
model of this sort is the so-called double-frequency sine- 
Gordon modelf£*& where an Ising transition was thought 
to be driven by two competing relevant perturbations. 
The description of the phase transition requires a non- 
perturbative scheme being able to identify correctly those 
degrees of freedom that remain massive and the low- 
energy ones that eventually become critical and undergo 
the transition. 

The same strategy was extended to a more compli- 
cated non- Abelian case — a two-leg spin- i Heisenberg an- 
tiferromagnetic ladder subject to a site-parity-breaking 
dimerization fieldi The non-trivial aspect of the non- 
perturbative approach adopted in this study of the 
strong-coupling limit of the model is that one has to 
preserve the SU(2) symmetry, which cannot be sponta- 
neously broken, in the Ising model language which appar- 
ently has only a discrete symmetry. This hidden symme- 
try was realized when using very non-local duality trans- 
formations for the coupled Ising models. An effective 
low-energy Hamiltonian depicting the ensuing quantum 
phase transition was derived and the transmutation of all 
physical fields at the infrared fixed point was identified. 
It was unambiguously shown that a quantum phase tran- 
sition occurs from the universality class of SU(2) level 2 
Wess-Zumino-Novikov-Witten (WZNW) model in the ul- 
traviolet limit to the universality class of SU(2) level 1 
WZNW model in the infrared limit. 



Another possible scenario was set up by introducing 
a staggered magnetic field, which explicitly breaks both 
SU(2) and bond-parity symmetries, to the generalized 
spin ladder^ Depending on which phase the model be- 
longs to in absence of the staggered field, either U(l) or 
Z 2 criticality was predicted under the effect of the stag- 
gered field. An interesting new phase that interpolates 
between the Haldane spin liquid phase and the sponta- 
neously dimerized phase was found beyond the transi- 
tions. The main features of this intermediate phase are: 
partially (transverse) coherent spin excitations and par- 
tially (longitudinal) non-vanishing string order parame- 
ter. 

As we can see, the physics of the generalized spin-i 
ladder under the effect of the dimerization or the stag- 
gered field is very rich. In this paper, we shall scrutinize 
these two models in a more general perspective and try 
to understand the phase diagrams as a whole. Our main 
attention is to be paid to the overall topology of the phase 
diagrams; while the physical properties in various phases 
are only to be mentioned briefly if available. 

The paper is organized as follows. In Sec. [HI we present 
a short overview of the continuous version of our models. 
In Sec. IHII we derive two sets of renormalization group 
(RG) equations and central charge formulas. These equa- 
tions are (numerically) solved and analyzed in Sec. II VI to 
bring about the phase diagrams of the continuous mod- 
els. The implication to the ladder systems is summarized 
and concluded in Sec. 



II. MODELS AND THEIR CONTINUOUS 
FIELD-THEORETIC MAPPINGS 

We consider a standard two- leg spin-i antiferromag- 
netic (J > 0) Heisenberg ladder (J±) generalized to in- 
clude a four-spin interaction (V)&i^ 



J ^""^ S a,n m S a,n+l + J± ■S'l/rt ' ^2,rt 



a— 1,2 n 
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superimposed by either a staggered magnetic field, model 
(I), 

Hh = hJ2 E(- 1 )" 5 a,„ , (2) 

a=l,2 n 

or a 7r-phase (relative) dimerization, model (II), 

a— 1.2 n 

Throughout the paper, we treat the models only in the 
weak-coupling regime, providing 

\J ± \,\V\,\h\{oi |A|)« J. (4) 

In this limit, the field-theoretic approach in continuum 
is of advance and power in describing the low-energy 
physics of the models. In fact, the resulting Hamiltoni- 
ans are represented by four interacting Majorana (real) 
fcrmions or, equivalent l y, fou r coupled Ising models. 

According to Refs.l^la l9llll we have the following map- 
pings in terms of Majoranas (£r,l) and Ising variables 

"7"^ gen — "^crit H - ^mass H~~ ^Ymarg •> (5) 

where the critical theory 

having an 0(4) symmetry is inherited from two decou- 
pled Heisenberg chains with SU(2)xSU(2) symmetry, 
v ~ JaU The mass terms result from part of the in- 
terchain coupling (Jj_) and the four-spin interaction (V), 
both of which have the same scaling dimension, d = 1: 

Wmass = - im s£,R^L ~ W™t£fl"£z, > ( 7 ) 

where the singlet (associated with £°) and triplet (£ = 
(C 1 !^ 2 !^ 3 )) masses are related to the parameters of the 
lattice spin model by 

m s = — 3cJj_ — c'V , m t — cJ± — c'V , (8) 

respectively. Here c and c' are some positive constants. 
The marginal terms are from the current-current part of 
the interchain coupling with scaling dimension d = 2: 

Wmarg = \ 9 l (Z R < L ) 2 + 92 .• £ L ) , (9) 

where c/i = | Jj_a — 7 and 32 = — \ J±o. — 7, with 7(> 0) 
the marginally irrelevant coupling in the chains. Ob- 
viously, the symmetry of the continuous model H ge n 
is 0(3) XZ2, corresponding to SU(2)xZ2 for the lattice 
model Q (Z2 arising from the interchange of two chains), 
unless J± = when the highest symmetry of 0(4) is 
reached, corresponding to SU(2)xSU(2). 



The field-theoretic counterparts of the staggered field 
PJl and the dimerization (j3J) are represented, respectively, 
by 

H h = (h/a)aia 2 fji3fjio , (10) 
Ha = (A/a)(Tia-2CT30-o , (11) 

where a with a magnitude order of the lattice constant 
is the short-distance cut-off of the theory. While Ha 
keeps the 0(3) invariance for the triplet sector intact, Hh 
breaks this symmetry down to 0(2) XZ2. This is a reflec- 
tion of the oriented field © lowers SU(2) to U(l)xZ 2 . 

In what follows, we shall address model (I): H gen _ + Hh, 
and model (II): H gcn +HA, in RG approach. Our purpose 
aims at drawing the phase diagrams by analyzing the 
infrared properties of the RG equations. 

III. DERIVATION OF THE RG EQUATIONS 

It is rather straightforward to establish the one-loop 
RG equations out of the operator product expansions 1 
(OPE's) between the perturbative operators^ For a crit- 
ical model perturbed by 

Hpcrt = — 7 KOi , (12) 
7ra z — ' 

i 

where Xi are dimensionless coupling constants of the cor- 
responding dimensionless operators Oi (with scaling di- 
mension d{), if Oi are normalized in such a way that 

(O i (z,z)O j (0M = SiAr^) , (13) 

(z = vr+ix being the time-space in complex coordinates) 
and the short-distance OPE's satisfy 

O t (z,z)O 3 (0,Q) nctjkf^LY ' k O k (0,0), (14) 
the one-loop RG equations for Xi are readily given by 
A fc = ^ = (2 - d k )X k - £ r./ X.X , + 0(X 3 ) . (15) 

y 

In addition, the central charge satisfying Zamolodchi- 
kov's decreasing theorem^ is perturbatively evaluated to 
this order to be 

C({X}) = C uv -3^(2-d fc )A 2 fe 

k 

+2]Tc J /'A 4 A J A fc + 0(A 4 ) , (16) 

ijk 

where Cuv is the central charge of the unperturbed crit- 
ical system (the ultraviolet fixed point). In the present 
case, Cuv = 2 is the central charge of two decoupled 
chains. 
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Under normalization JT3J, the OPE coefficients Cy k are 
related to the three-point function (OiOjOk), and thus 
are complete symmetric among the indices (ijk)*^ Not 
only this procedure simplifies the calculations, but also 
is it essential to give rise to a correct form for the central 
charge ifTHjl . 



A. RG equations for model (I) 



Since for this model, 

' '-pert ' L n 



7~i-h 



(17) 



and noticing TLh breaks the 0(3) symmetry in the triplet 
sector, we have to divide the triplet mass m t and the 
marginal coupling constants g\ and <?2 into doublets and 
singlets. Denoting the mass bilinear (or energy density) 
of the Ising model, — i^^, we write down: 



H 



pit = ~ m t( e i + £ 2) - m s t s 3 - m s e 

-fli£i£2 -5i(ei + ^2)^3 
-giiei + e 2 )e Q - g s 2 e 3 e Q 



v \ - 



na 



2 



(18) 



where the dimensionless operators in both fermion and 
Abelian boson representations ("+/—" associated with 
sectors (l,2)/(3,0), respectively): 

0\ = -V2na(ei +e 2 ) = -V2 cos V4n$ + , 

2 = -2nae 3 = - cos \/47r$_ + cos V47rO_ , 

03 = — 2-iraEQ = — cos V47t < I > - — cos V^t©- ; 
C 4 = -(27ra) 2 £i£ 2 = 47ra 2 <9$ + <9$+ , 

5 = -2V2(7ra) 2 ( £l +e 2 )e 3 
= -v2cos V^7r ( I > +( cos \/47r$_ — cos \/47r6-) , 

6 = -2V2(7ra) 2 ( £l +e 2 )eo 

= -V2 COS \/ 4tT$_)_ f COS \J 47T < £ , _ -+- COS \/ 47T@_) , 

C 7 = -(27ra) 2 e 3 eo = 47ra 2 <9$_<9$_ ; 
Og = 2cricr 2 /Z3^o = 2 sin \A" ( I > + cos y^^- • (19) 

The scaling dimensions, ^1,2,3 = 1, (^4,5,6,7 = 2, and d s = 
h. Accordingly, the dimensionless coupling constants: 
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It is easy to calculate all the OPE coefficients by using 
either fermion or boson representation. For example, we 

O s (z, z)Os ~ + ^f x - \{0 2 + 3 ) 



(21) 



Therefore, we read off 
1 



C88 



>/2' 



C88 5 — C88 



2V2 ' 



, 7 = -j. (22) 



The algebra in this way is closed, and we list all the 
coefficients in Table [I] 



TABLE I: OPE coefficients cy fc ( given in the parentheses, in 
front are indices k). 



i\j 1 2 3 4 5 6 7 8 



4(-l) 5(-l) 6(-l) 1(-1) 2(-l) 3(-l) 



1(-1) 



3(-l) 8(4) 
1(-1) 2(-l) 8(-|) 



5(-l) 7(-l) 

6(-l) 7(-l) 

1(-1) 5(-l) 6(-l) 8(4.) 

2(-l) 1(-1) 5(-l) 4(-l) 7(-l) 6(-l) 8(4^) 

3(-l) 1(-1) 6(-l) 7(-l) 4(-l) 5(-l) 8(4^) 

3(-l) 2(-l) 6(-l) 5(-l) 8(4) 

8(73) 8(4) 8(4) 8(-|) 8(4.) 8(4.) 8(4J ( C 88 fe )° 



8 C8 s fc = l(^), 2(-i), 3(-|), 4(-I), 5(4^), 6(4^), 7(-$). 

Now, it follows from Eq. (|15|) that the RG equations 
for the couplings 77 (instead of A): 

1 , s 1 2 

% = r h + + + mm) - jVs > 

N 1 2 

772 = m + ^(2mv5 + mvi) + t% » 
% = % + -jymm + 772777) + t% ; 
% = H + \(vi + vl) + \vl , 

!/ N 1 2 

775 = 2^17/2 + -(774775 + 776777) - ^77 8 , 

1, , 1 o 

776 = 2771773 + -(774775 + 775777) - -77 8 , 

1 2 

777 = 2773773 + 775776 + -77 8 ; 

3 1. 
t?8 = 2 7? s ~ 7 ?i^8 + + 773)778 

-4(775 + 776)778 + 7^(774 + 777)778 . (23) 

(20) 

The perturbative central charge (|16|) for the present 
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model in terms of 77 is 
3 



and the central charge to 



C(M) = 2-|(2 77l 2 + % 2 + % 2 + |% 2 ) 

3/ 11 



3/„ 1 
"8 V m ~ m ~ V3 ~ 4 
1 1 1 
+ 2^5 + 2% - ^J^g 



(24) 



B. RG equations for model (II) 

Comparing with the previous model, the only differ- 
ence now is replacing Hh by Ha, i-e., 



08 = 2CT1CT2C3C0 = 2 sin v / 7T c I > + sin v / 7r ( E > - 
% _ 7rA 
2 ~~ 2t; ' 



A; 



(25) 
(26) 



and Eq. (|2"T]l is modified to be 



O s (z,z)O s 



a 1 

i 



0i 



t 



2V2 



(o 5 



2 

6 ) 



(o 2 + o 3 ) 
1 



(04 + 7 ) 



yielding 



4(-l)>5(-^),6( 



2^2) 



7 



(27) 



(28) 



However, as we have known, the symmetry of model 
(II) is higher than that of model (I) — the 0(3) symme- 
try in the triplet sector [i.e., the SU(2) symmetry in the 
original lattice model] should remain. In other words, 
the division of the triplet into doublet and singlet is fu- 
tile in the present case. This fact is reflected in the RG 
equations also. From Eqs. 1)23(1 . by reversing the signs of 
the coefficients involving ess 2 , ess 3 , ess 5 , and ess 6 , three 
of the eight equations are redundant, implying 771 =772, 
Vi = V5, and = V7- As a result, the RG equations 
degenerate to 



m 


= 771 + - (2771774 


^ 1 2 

+ mv6) - 4% 


m 


3 

= m + - 


1 2 


m 


= 277? + + 




m 


= 2771773 + 774776 


1 2 


m 


= 2^ ~2 {3m 





(^4 + ?76)778, (29) 



C(W) 



% 2 



3 2 \ 
2 % J 



3 
"16' 

-§(* 



% 2 



(30) 



IV. SOLUTIONS TO THE RG EQUATIONS 
AND THE EMERGING PHASE DIAGRAMS 

To analyze the infrared behaviors of our present mod- 
els, in this section, we numerically integrate these two 
sets of RG Eqs. (|2"3|) and (|2"9"|) to depict generic phase 
diagrams for the spin ladder under a staggered field or 
dimerization. On account of the fact that for a system in 
a massive regime, the role of marginal coupling is usually 
exhausted by renormalizing the mass and velocity, it is 
reasonable to take no thought of these marginal perturba- 
tions at the outset and argue that the overall qualitative 
feature remains correct although the position and shape 
of phase boundary may not be exact. Let us begin with 
a known system — the generalized spin ladder (|T)l^ — as an 
example. A variant of this model includes four-spin ring 
exchange, which is believed to be relevant to some exotic 
properties in real ladder materials £L1L1!L 

In the following we denote the bare ultraviolet cou- 
plings by 77'°) = 77(L = 1) and the renormalized infrared 
couplings by t^ 00 ) = r\(L = 00). 



A. Phase diagram of the generalized spin ladder 

In this case, the RG equations are trivially 771 = 771 and 
% = V3- Unless the initial value (bare coupling) vanishes, 
77^ = or 773°'' = 0, the system is always renormalized to 
some strong-coupling massive phase as L — > 00. Because 
of relations in Eqs. 77^ = ^ and t?^ = ^f-, 

we conclude that to* = and ra s = are phase transi- 
tion lines on which the model becomes critical. They be- 
long to the universality classes of critical SU(2) 2 WZNW 
model and Z 2 Ising model, respectively^ As a base for 
our further exploration, we draw the phase diagram of 
the generalized spin ladder Q in Fig. ^ 

According to Ref. there are four massive phases di- 
vided by two critical lines. The second and fourth quad- 
rants of the figure are two Haldane (H^) spin liquid 
phases, and the first and third quadrants belong to spon- 
taneously dimerized (sDw) phases. The "+/— " sectors 
are actually related by duality of reversing the signs of the 
masses. The axis of abscissas represents a critical model 
with an SU(2) 2 WZNW universality class (L:SU(2) 2 ), 
and the axis of ordinates a Z 2 (Ising) one (L:Z 2 ). Al- 
though the Haldane phases and the dimer phases share 
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V 



L: Z, 



S: H (+) 
I I I I I I 



S: sD 



(-) 



S: sD 



<+) 



L: SU(2) 2 



r| ( 3 0) ~ m s 



S:H 



(-) 



the central charge. For this simple model we have 



(31) 



Considering this is a perturbative one-loop result and 
should not be taken seriously when 77 approaches unity, 
we assume that we can mimic the true fixed point prop- 
erty by introducing a cut-off in 77. If we are allowed to 

choose this cut-off to be «/|, the correct central charges 

for both SU(2)2 and Z2 can be recovered, since for the 

I , and for the lat- 



former when (771,773) = (0, y|), C* 



ter when (771,773) = (y|,0), C* = |. The existence of 

the SU(2)2 criticality was recently confirmed by a direct 
estimate of the numerical value of the central charged 



FIG. 1: Phase diagram of the generalized spin ladder. There 
exist two critical lines L: SU(2)2 and L: Z2, and four massive 
sheets S: H (±) and S: sD (±) . Notice the dimers in the Haldane 
phases are only symbolic rather than a true order. 



B. Phase diagram of model (I) 

Model (I) is the generalized spin ladder in a stag- 
gered field (JSJ . The RG equations regarding the relevant 
perturbations now become 



the feature of non- vanishing topological string order pa- 
rameter in common and thermodynamically they are in- 
distinguishable, the symmetry of the ground state and 
the behavior of the dynamical susceptibility arc com- 
pletely different. 

(i) S: HW phase = +00, rj^ = -00]: Neither 
parity (site or bond) symmetry nor translational symme- 
try is broken. There is a coherent peak for the dynamical 
spin susceptibility at<7_L = 7r,<7~7r,w^ \mt\, charac- 
terizing coherent S = 1 single-magnon excitations. 

(ii) S: H<-> phase = -00, rj^ = +00]: No 
symmetry breaking (except for the hidden string order). 
The peak is at q± = 0, q ~ n, lu ~ \m t \ instead. 

(iii) S: sD(+) phase [rj^ = +00, = +00]: 
The discrete site-parity and translational symmetries are 
spontaneously broken. The corresponding order parame- 
ter is ({-l) n {S hn -S 1 , n+1 + S 2 , n -S 2 , n+1 )) ± (0-phase). 
No coherent peak in the spin susceptibility (square-root 
singularity only). 

(iv) S: sD^") phase [rj^ = —00, = —00]: The 
discrete symmetries are broken with an order parameter 

((-l) n (Sl 1 n-Sl,n+l - S 2 ,n-S 2 ,n+l)) ¥= (^-phase). No 
coherent peak either. 

We also show in the figure the directions of the lattice 
model parameters J± and V via relation ©, although 
the exact scales for the coordinates are not known. We 
find that the standard ladder with positive, antiferro- 
(negative, ferro-) J± is in Haldane H(+) (HH) phase. 
The pure four-spin coupling model with J± = is in the 
sD(~) phase when V > or in the sD' + ) phase when 
V < 0. This is a special case with a higher symmetry. 

It is also interesting to check in RG sense the flow of 





= Vi - 


1 2 

1% ' 
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= m + 


1 2 

1% ' 




m 


= V3 + 


1 2 




m 


3 

= - 2 m 


- mm h 





(32) 



Evidently, the equations are invariant under the trans- 
formation 77 8 — > —77g, which corresponds to h — > —h in 
the original spin model. We need only deal with the 
positive case. At first glance the system has two mathe- 
matical fixed points: (a) (77* , 772, 773, 77I) = (0,0,0,0) and 
(b) (|, — |, — |, y/S). The one at the origin is an unstable 
fixed point recounting two decoupled chains as our start- 
ing point; while the fixed point (b) is fallacious. This is 
due to the fact that 771 and 772 are not completely inde- 
pendent, i.e., initially 77^ = r/^ = 2HS ) rendering the 
fixed point (b) physically unaccessible. However, the val- 
ues in (b) do have full play in the the structure of the 
phase diagram. 

The emerging phase diagram of model (I) is shown in 
Fig. 121 It is notable that all the interesting phases lie 
in the first three quadrants of the base and the fourth 
quadrant is in some sense rather barren. This is con- 
trolled by the "fixed point" (b) mentioned above. To 
make the structure clearer, we anatomically plot figures 
in two intersecting planes (see Fig. Various phases 
are labeled in the figures, where "B" denotes a block or 
bulk space that is always massive, "S" stands for a sur- 
face or sheet as phase boundary separating two different 
phases, which can be critical or first-ordered, and "L" is 
an intersectional line of two sheets. 
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l 0.1 




Tl<°Wl< 2 0) ~ m t 



where ij^'s are bare-couplings as the initial values of the 
differential equations rj(°> = rj(L = 1 ). 

Now the system is subjected to a staggered field. It is 
found when the field is increased to 



0.0354863 , 



(34) 



the system undergoes a second-order phase transition 
with a U(l) universality class (or Gaussian fixed point). 
The behaviors of the RG flows in the vicinity of this 
transition is shown in Figs. 0] (a), (b), and (c). Below 



FIG. 2: Phase diagram of the generalized spin ladder in a 
staggered field. The massive phases occupying different blocks 
are separated by interfaces either second-ordered (enclosed by 
solid line) or with the transition nature unknown from the RG 
equations alone (dashed line). 



0.15 



0.1 



0.1 



Si CO 



0.05 



(a) r,(°)=-0.2 



B:H H 




B:H <+) 



S:Z .v- I /, ,S:U(1) 
B-sD H L:SU ' 2 '- 



-0.2 



-0.1 



0.1 



0.2 



r,(o)=r|(°> 



0.15 



0.1 



0.05 





(b) r|< >=r|< )=0.2 
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FIG. 3: Sectional plot of the phase diagram of the generalized 
spin ladder in a staggered field: (a) ^g ' is a constant as —0.2. 
(b) j^ ' = r}^ is a constant as 0.2. 



Now let us interpret in some detail the obtained phase 
diagram. We take the example with a parameter range 
where two consecutive phase transitions occur when in- 
creasing the staggered field — one belongs to the U(l) uni- 
versality class; the other's property is not determined by 
the RG theory itself (probably first-ordered). Suppose a 
generalized spin ladder offers such a kind of parameters: 



7* 0) = Vf ] = 
< = -0.2 , 



0.01 



(33) 
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(f) . 
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FIG. 4: Parametric plot of the RG flows across two phase 
transitions. The figures in the left column describe the flows 
in the vicinity of the U(l) transition, where the solid lines 



correspond to ?7g 



(0) 



(oh 



the dashed lines to ?7g 



(()) 



(% )d -0.001, and the dotted lines to ' = (r/s )a +0.001. 
The figures in the right column are in the vicinity of the sec- 
ond (unidentified) transition, where the solid lines correspond 
to 7^, 0) = (r^ 0) ) c2 , the dotted lines to rj ( s 0) = (r7 8 0) ) C 2 - 0.001, 
and the dash-dotted lines to ?7g ' = (^g '^ + 0.001. Notice in 
panel (f) the dash-dotted line for r/3(L) eventually approaches 
+oo after an upturn that is not shown in the figure. 

the U(l) transition the system is in the bulk massive 
phase B: H( + ) characterized by the couplings renor- 
malized to (r)[°°\ rj^ > V3°°\ vi^) — (+oo, +oo, — oo, 0) 
(the dashed lines in the figures). At the transition 



the infrared property becomes (r]i°°\ "q^^ , 773° 



(0, +00, — 00, 0) (the solid lines in the figures). Since 
771 is proportional to the split doublet mass, mf, it is 
bound to attribute this transition to a U(l) character. 
Now we are in the bulk phase of B: I characterized by 
(—00, +00, —00, 0) (the dotted lines in the figures). This 
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is an intermediate phase interpolating between the Hal- 
dane phase and the spontaneously dimered phase. The 
most distinctive features of this phase are that it has 
only transverse coherent spin excitations and longitudi- 
nal string order parameter. (For detailed properties, cf. 
Ref.H) 

If the staggered field is further increased to reach the 
point 

h c2 ~ (ri { s ] )c2 = 0.0560487, (35) 

the second transition occurs. The RG flows around 
this transition is shown in Figs. 0] (d), (e), and (f). 
At the transition the couplings are renormalized to 
(— oo, +oo, — oo, y/S) (the solid lines in the figures). After 
the transition the system enters the B: H^~J phase char- 
acterized by (— oo, +oo, +oo, oo) (the dash-dotted lines 
in the figures). Unfortunately, we can tell very little 
from the present scheme regarding the property of this 
transition. It may be a first-order transition or even a 
crossover. 

In this way, we are able to classify most of the phases 
and transitions by the infrared properties inferred from 
the RG equations. We summarize the main results as 
follows. 

a. Bulks : 





(oo) 

m 


(oo) 


(oo) 

V 3 


(oo) 

Vs 


B:sDW 


+oo 


+oo 


+oo 





B:sD("> 


— oo 


— oo 


— oo 





B:H(+) 


+oo 


+oo 


— oo 





B:H<~) 


— oo 


+oo 


+oo 


00 


B:I 


— oo 


+oo 


— oo 






(36) 



Comparing with the zero-field case (previous subsection), 
we find that three massive phases sD^ and H^ + ) are es- 
sentially passed from their h = counterparts directly 
except for the broken symmetry rji + 1 r/2- However, 
B: H^ - ) phase is different from S: H^ - ) in the sense that 
772 ~~ > instead of — oo and 1]%°°^ — > oo instead of 
0. Obviously, B: I is a totally new phase which has no 
zero-field counterpart. 
b. Sheets : 





(oo) 

m 


(oo) 

m 


(oo) 


(oo) 

% 


S:U(1) [HW-I] 





+oo 


— oo 





S:Z 2 [sD(")-I] 


— oo 





— oo 





S:Z 2 [H(+)-sD(+)] 


+oo 


+oo 








S:Z 2 [sD(-)-H(-)j 


— oo 


— oo 


— oo 


V3 


S: [I-H(-)] 


— oo 


+O0 


— oo 


V3 


S: [H(+)-H(-)] 


+oo 


+O0 


— oo 


V3 


S: [sDW-H(-)] 


+oo 


+00 


+00 


V3 



(37) 



It is perspicuous that the critical surfaces split from the 
SU(2)2 critical line and having U(l) and Z 2 criticalities 
correspond to the infrared fixed points with (771 , 77s) oc 
(mf,h) — > (0,0) and (772,77s) oc (m*,h) — ► (0,0), respec- 
tively. The other Z 2 criticality is trivially same as and 



connected to the one when h — 0, (773,778) ^ ( m s,h) — > 
(0,0). At least one of the four unknown boundaries 
(dashed lines in Fig. and 0J) can be identified as Z 2 , 
i.e., S: [sD( - )-H' - )]. The reason is to be shown soon, 
c. Lines : 





(00) (00) 


(00) 
m 


(00) 

Va 


L:U(1) [Hl+I-I-HH] 


+| +00 


—00 


V3 


L:Z 2 [sD(-)-I-H( - )] 


-00 -§ 


—00 


V3 


L:Z 2 [H<+)-sD(+)-H(-)] 


+00 +00 


3 

4 


V3 



(38) 



These lines are intersections of two surfaces. The special 
values here are associated with the "fixed point" (b) of 
Eqs. @2J. 




V a 



FIG. 5: A top- view plot of the phase diagram of the general- 
ized spin ladder in a staggered field. The true phase bound- 
aries in the base (h = 0) should be extended as shown by the 
arrows. 

To close our discussion, we show in Fig.[S]a top- view of 
the phase diagram. The U(l) zone is confined in a wedged 
area demarcated by the SU(2) 2 line 77^ = rj^ = 

and rii = m — ~^ r l^ > 0, i.e., m t = —\ m s > 0. 
The latter condition amounts to V = and J± > 
via relation (JHJ), coinciding with the location of a stan- 
dard spin ladder with antiferromagnetic interchain cou- 
pling. The adjacent Z 2 region is delimited by condition 
rj J ' = 772°'' = hvs < 0. The trivial Z 2 region is within 
the Z 2 line 77^°' = and 77^ = t/^ ' = -3^ 0) > 0. The 

unidentified sheet have 77^°'' = = 773 as its bound- 
ary in the base. This is a model with J± = 0, i.e., only 
four-spin interchain coupling. However, one might have 
noticed that there is something unnatural here. Consid- 
ering the internal relation between the B: H^ - ) phase in 
the bulk and the S: H^ - ) phase in the base, we suggest in 
the true situation the boundary of this sheet should bend 
over to connect to the critical lines in the base (as shown 
by the arrows in Fig. |SJ) ■ By combining the "integrating 
out" method and the study of the bosonized model, we 
observe that a Z 2 surface, connecting two axes, covers 
the whole third quadrant in Fig. [3] 
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C. Phase diagram of model (II) 



Now we turn to Model (II) , the generalized spin ladder 
JTJ under dimerization J3}. It follows from Eqs. I|29[l by 
retaining only the relevant couplings: 



SU(2) 



m 


= Vi - 


4 






1 


>h 


= m- 


4 




3 




is 


= - 2 m 





(39) 



This is a model with an unbroken 0(3) symmetry. As 
before, we have an invariance rjs — > — r/g (i.e., A — > —A). 
In contrast with the previous model, now the system has 
really a fixed point (77* , 773, r^g) = (1,1, v3) besides the 
trivial one at the origin. This fixed point occurs at 771 = 
773, or equivalently to* = m s , corresponding to the pure 
four-spin coupling model ( J±_ — 0) when the system has 
a higher SU(2)xSU(2)«0(4) symmetry. This is also the 
condition for the line where the SU(2)i sheet merges the 
Z2 sheet into the SU(2)2 sheet. The RG flows around 
this fix point are shown in Fig. where we choose the 




RG flows to the fixed point (77* , 773 , 7?|) = (|, |, v3) 



FIG. 

(solid line) at r][ 0) = r?f ] = 0.1 and (rj ( ° ] 



below (^ 0) = (77^) 



0.137410, and 
0.001, dashed line) and above (T/g ' = 



(77g°') c + 0.001, dotted line) the transition. 



initial values for the base quantities: r}\ — =0.1. 
The initial value of the dimerization for the fixed point 
is (^ 0) ) c = 0.137410. 

The phase diagram of this model is plotted in Fig. 
and Fig. is for the sectional pictures of the same phase 
diagram, where we have labeled various phases. As be- 
fore, the classification of phases and criticalities is asso- 
ciated with the infrared properties of the system resulted 
from the RG equations. 




r|< 0) ~m s 0.1 



0.1 r)', ' ~m t 



FIG. 7: Phase diagram of the generalized spin ladder under 
7r-phase dimerization. 
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B:sD <+) 


B:H (+) 


L:Z 2 { 


^ S:Z 2 



-0.2 



-0.1 



0.1 



0.2 



TIP) 



0.4 



0.3 



0.2 



0.1 



(b) r|<°)=0.2 




B:D ( -) 

S: 




B:II ^ 

^SC'jy L:SU(2) 2 


B:sD (+) 


- S:SU(2) 2 



-0.2 



-0.1 



0.1 



0.2 



T1V 



(0) 



FIG. 8: Sectional plot of the phase diagram of the generalized 
spin ladder under dimerization. 



a. Bulks 





(00) 

m 


(00) 

V3 


(00) 


B:H(+) 


+ 00 


— 00 





B:H(") 


—00 


-|-oo 





B:sD(+) 


+00 


+CO 





B:D(") 


—00 


— 00 


OO 



(40) 



Three massive phases H^) and sD(+) are same as their 
A = counterparts. Since the symmetry of phase sD^ + ) 
is different from the explicit dimerization perturbation 
©, it is still spontaneously dimerized. However, B: D'~' 
phase is different from S: sD' _ \ because now the system 
is explicitly dimerized with the same symmetry of the 
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perturbed Hamiltonian. In the infrared limit T]g — > oo 
instead of 0. It is also noticed from the figures that phase 
B: H^ - ) is restricted to a rather small region in the phase 
space. This may not be the real situation (see discussions 
at the end of this section). 
b. Sheets : 





(oo) 

Vi 


(oo) 

m 


(oo) 

Vs 


S:SU(2) 2 [aDW-Hl-J] 





+O0 





S:Z 2 [gDW-HW] 


+oo 








S:SU(2)i pW-D(-)] 


+oo 


— oo 


V3 


S: p<->-D<->] 


— oo 


+oo 


V3 



Like the case before, the identification of the SU(2)2 and 
Z 2 criticalities is directly due to the fact that in addi- 
tion to the dimerization (1]$°°^), the triplet (f][°°^) and 

singlet (7J3 ) masses are renormalized to vanish, respec- 
tively. However, the reference of the interface between 
phases H^ + ) and D^ - ) to an SU(2)i criticality is a bit 
subtle. This conclusion is drawn from other studies of the 
standard ladder under dimerization£i2i22i2i In Ref. HH 
a standard spin-^ ladder with ferromagnetic interchain 
coupling under 0-phase dimerization (the dual case of 
our present model) was semiclassically mapped onto an 
0(3) nonlinear sigma model with some topological an- 
gle. The SU(2)i criticality is realized when this angle 
equals 7r. On the other hand, an effective Hamiltonian 
was explicitly derived in Ref. []] which describes the same 
transition in terms of a single spin-i chain with bond 
alternation, and the transmutation of all physical quan- 
tities from the ultraviolet fixed point to the infrared fixed 
point was established. We have only one unidentified sur- 
face, S: [Hy-)-D(-)], in this model, 
c. Line : 





( OO ) ( OO ) ( OO ) 

m m vs 


L:SU(2) 2 [sD(+)-H(+)-D(-)-H(-)] 


33 /0 

4 4 V<5 



(42) 



The only critical line in this model represents the fixed 
point of Eqs. 1)39(1 . the RG flow for the point on which is 
shown in Fig. This is a very special line with a higher 
0(4) symmetry: it is not only the ridge line where all 
sheets meet, but also the common boundary shared by 
all four phases (blocks) of the model. 

Since the obtained phase diagram is grounded on the 
relevant part of the one-loop RG equations, the posi- 
tion of the phase boundaries may not be so correct. 
For instance, considering the intimacy between the bulk 
B: D^~' phase and the S: sD' - ) phase in the base, most 
plausibly the SU(2)i sheet is connected to the SU(2) 2 line 
of A = 0, and the boundary for the unidentified sheet is 
pushed to connect to the Z 2 line in the base (see Fig- - 
So the space for the bulk B: H^ - ) phase is considerably 
enlarged. This is supported by the fact that, as we have 
mentioned, the (Jj_ > 0) standard ladder (V = 0) under- 
goes the SU(2)i transition under the (7r-phase) dimeriza- 
tion. So the SU(2)i sheet cannot terminate at the loca- 




V A 

FIG. 9: Top- view plot of the phase diagram of the generalized 
spin ladder under dimerization. The true SU(2)i sheet should 
connect to the SU(2)2 line in the base, and the unidentified 
sheet to the Z2 line in the base (as shown by the arrows). 

tion which is a result from the approximate RG Eqs. I|39|) 
as shown in Fig. [7| 

V. SUMMARY AND DISCUSSIONS 

Based on the one-loop RG equations and the assump- 
tion that the relevant perturbations govern the relative 
structure of various phases in the parameter space, we 
have speculated from the infrared properties of these 
equations the phase diagrams of two related spin ladder 
models in the weak-coupling regime. 

In model (I), the external staggered field breaks the 
SU(2) symmetry down to U(l) and Z 2 , so that the criti- 
calities at most have the same symmetries and no room 
for SU(2) to live. Indeed, we find both U(l) and Z 2 criti- 
calities come forth with the introduction of the staggered 
field (see Fig. and |5j) , which actually evolve from the 
underlying SU(2) 2 criticality in absence of the field. The 
present result confirms the picture we obtained through a 
somewhat different approach in Ref. |8j In that study, we 
literally integrated out the fast degrees of freedom associ- 
ated with the singlet modes close to the SU(2) 2 line and 
derived an effective action which describes the triplet sec- 
tor with anisotropy induced by the staggered field. The 
vanishing of the split effective triplet mass manifests the 
criticality accordingly. The new observation is that the 
intermediate phase (B: I), which interpolates between the 
Haldane phase (B: H^ + ') and the dimer phase (B: sD(~)), 
is restricted to a very limited region, and there should ex- 
ist another (unidentified) transition which separates B: I 
phase from the large-/i (Haldane) phase (B: H^ - )) which 
occupies most of the space in the phase diagram. From 
the infrared property we see that B: I phase is quite 
unique entirely due to the staggered field, and there is no 
resemblance preexisting when h = 0. 

Comparing with model (I) , model (II) is relatively sim- 
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pie in the sense that the SU(2) symmetry is unbroken and 
the 0(3) symmetry in the triplet sector remains, and 
the number of RG equations is thus reduced. All four 
bulk phases can basically be traced back to their coun- 
terparts when A = 0. None the less, as we have seen, the 
phase diagram of this model is by no means uneventful 
(see Fig. \J\ and |5J). The phenomenon extremely intrigu- 
ing is the merging of SU(2)i and Z 2 critical sheets into 
SU(2) 2 sheet. The information encoded in the RG equa- 
tions reveals that this occurs at some critical dimeriza- 
tion strength A C (V) when the system has a higher 0(4) 
symmetry or SU(2)xSU(2) in the original lattice model, 
i.e., Jj_ = 0. The merging line with SU(2)2 criticality is 
shared by all four massive phases of the model. 

The benevolent aspect of our present theory is that 
we are able to assort different phases according to their 
infrared properties from the RG flows, and the identifi- 
cation of the critical sheets that emerge from the under- 
lying critical lines in the base should be reliable. These 
critical sheets, including two Z 2 and U(l) in model (I) 
and SU(2) 2 and Z 2 in model (II), can also be verified 
to exisli^ by means of "integrating out" procedure em- 
ployed in Ref. 8, although the terminating positions of 
these sheets are unknown by this approach. However, 
as one may have perceived, the insolvency of the present 



theory is that not all of the transitions can be attributed. 
Even worse is the phase boundaries resolved from the 
one-loop equations may not be accurate. (See our sug- 
gestions for the modification of the phase diagrams in 
Fig. [5] and |5|) Therefore, although the topology of the 
phase diagrams we believe has been overall and qualita- 
tively captured, the unwavering pin-pointing of the phase 
diagrams has to resort to other, say numerical, methods. 
For example, it is predicted here for model (II) of pure 
four-spin interaction (V), when increasing the dimeriza- 
tion, there is only one second-order (SU(2) 2 ) transition 
which connects directly two dimerized phases with mu- 
tually mismatched orders. If adding a bit positive J± 
coupling, one would expect two consecutive transitions: 
Z 2 and SU(2)i in a row. 



Acknowledgments 

It is a pleasure to acknowledge the helpful discussions 
with P. Fulde, K. Penc, R. Narayanan, T. Vekua, and 
M. Nakamura. Special thanks are due to A. Nersesyan 
for calling the problem to my attention. This work was 
supported under the visitors program of MPI-PKS. 



1 See, for a concise review on conformal field theory, P. 
Ginsparg, in Fields, Strings and Critical Phenomena, ed. 
by E. Brezin and J. Zinn- Justin (Elsevier Science, Amster- 
dam, 1989). 

2 See, for a review on spin chains, I. Affleck, in Ref. [l| J. 
Phys.: Condens. Matter 1, 3047 (1989). 

3 See, for a recent book on bosonization, A.O. Gogolin, A. A. 
Nersesyan, and A.M. Tsvelik, Bosonization and Strongly 
Correlated Systems (Cambridge University Press, Cam- 
bridge, England, 1999). 

4 See, e.g., M.A. Continentino, Quantum Scaling in Many 
Body Systems (World Scientific, Singapore, 2001). 

5 G. Delfino and G. Mussardo, Nucl. Phys. B 516, 675 
(1998). 

6 M. Fabrizio, A.O. Gogolin, and A. A. Nersesyan, Phys. Rev. 
Lett. 83, 2014 (1999); Nucl. Phys. B 580, 647 (2000). 

7 Y.-J. Wang and A. A. Nersesyan, Nucl. Phys. B 583, 671 
(2000). 

8 Y.-J. Wang, F.H.L. Essler, M. Fabrizio, and A. A. Ners- 
esyan, Phys. Rev. B 66, 024412 (2002). 

9 A. A. Nersesyan and A.M. Tsvelik, Phys. Rev. Lett. 78, 
3939 (1997). 

A.K. Kolezhuk and H.-J. Mikeska, Phys. Rev. Lett. 80, 
2709 (1998). 

1 D.G Shelton, A. A. Nersesyan, and A.M. Tsvelik, Phys. 
Rev. B 53, 8521 (1996). 

2 Rigorously speaking, the interchain coupling J± shifts the 



triplet and singlet velocities differently. However, we ne- 
glect this difference under the condition (J1J. 

13 See, e.g., J.L. Cardy, in Ref. H A.W.W. Ludwig and J.L. 
Cardy, Nucl. Phys. B 285, 687 (1987). 

14 A.B. Zamolodchikov, JETP Lett. 43, 730 (1986) [Pis'ma 
Zh. Eksp. Teor. Fiz. 43, 565 (1986)]. 

The definitions and conventions we take may be found, 
e.g., in Appendix A of Ref. H 

16 S. Brehmer, H.-J. Mikeska, M. Muller, N. Nagaosa, and S. 
Uchida, Phys. Rev. B 60, 329 (1999); M. Matsuda, K. Kat- 
sumata, R.S. Eccleston, S. Brehmer, and H.-J. Mikeska, 
ibid. 62, 8903 (2000); T.S. Nunner, P. Brune, T. Kopp, M. 
Windt, and M. Griininger, ibid. 66, 180404(R) (2002); 

17 K. Hijii and K. Nomura, Phys. Rev. B 65, 104413 (2002). 

18 M. Muller, T. Vekua, H.-J. Mikeska, Phys. Rev. B 66, 
134423 (2002); T. Hikihara, T. Momoi, and X. Hu, Phys. 
Rev. Lett. 90, 087204 (2003); A. Lauchli, G. Schmid, and 
M. Troyer, Phys. Rev. B 67, 100409(R) (2003). 

19 K. Totsuka and M. Suzuki, J. Phys.: Condens. Matter 7, 
6079 (1995). 

20 M.A. Martin-Delgado, R. Shankar, and G. Sierra, Phys. 
Rev. Lett. 77, 3443 (1996); M.A. Martin-Delgado, J. 
Dukelsky, and G. Sierra, Phys. Lett. 250A, 430 (1998). 

21 D.C. Cabra and M.D. Grynberg, Phys. Rev. Lett. 82, 1768 
(1999). 

22 A. A. Nersesyan and Y.-J. Wang, unpublished. 



